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We explore multiband effects on d-wave superconducting instabilities driven by electron-electron 
interactions. Our models on the two-dimensional square lattice consist of a main band with an 
extended Fermi surface and predominant weight from d^2_y2 orbitals, whose orbital character is 
influenced by the admixture of other energetically neighbored orbitals. Using a functional renor- 
malization group description of the superconducting instabilities of the system and different levels 
of approximations, we study how the energy scale for pairing and hence the critical temperature 
is affected by the band structure. We find that a reduction of orbital admixture as a function of 
the orbital energies can cause a Tc enhancement although the Fermi surface becomes more curved 
and hence less favorable for antiferromagnetic spin fluctuations. While our study does not allow a 
quantitative understanding of the Tc differences in realistic high-Tc cuprate systems, it may reveal 
an underlying mechanism contributing to the actual material trends. 

PACS numbers: 74.20.-z, 74.72.-h, 74.62. Bf 



I. INTRODUCTION 

With the advent of the new iron superconductors,"'^ 
there is now a second class of superconducting materials 
with high Tc's besides the layered cuprates. As opposed 
to the cuprates, the iron pnictide superconductors have 
more than one electron band at the Fermi level, and the 
present theoretical works to understand these materials 
clearly indicate that this fact makes the systems far more 
complex. In particular, theorists have undertaken efforts 
to relate the observed differences in the transition tem- 
peratures and in superconducting gap structures through 
the families of iron superconductors to details in the crys- 
tal and hence electronic structurei^"— Similar ideas were 
then applied to the cuprates as welli^ These works tackle 
the full complexity of many-orbital problem with approx- 
imate many-body techniques. The results are interest- 
ing and promising, as they show the existence of various 
tuning parameters for Tc- On the other hand, we feel 
that the problem of correlation-driven superconductivity 
should also be approached from a constructive point of 
view using simpler models, by asking what changes occur 
if one takes a separated Fermi surface with a supercon- 
ducting instability, and adds to it the orbital character of 
the band, or allows other bands to come close in energy. 
The main question is whether one can find explainable 
trends that can be used as guide lines in a search for 
higher Tc's or other desired properties. This is one of the 
motivations for the study described below. 

Another direct motivation for considering electron- 
electron-interaction-driven pairing in multiband models 
is the Tc trend in the high-Tc cuprates. As pointed 
out ten years ago by O.K. Andersen's group there ap- 
pears to be a positive correlation between the experi- 
mental Tc's and the theoretically derived second-nearest- 
neighbor hopping parameter t' or a related parameter r. 
Here, the higher Tc's occur for a rounder Fermi surface. 



i.e., for larger t' and larger r. This contradicts at least 
the naive expectations in a spin-fluctuation-induced pair- 
ing scenario, where a smaller amount of nesting leads to 
a weaker pairing interactions, and hence a round Fermi 
surface with larger t' should, at least over some parameter 
range, have a smaller Tc. This thinking is purely based on 
the geometry of the Fermi surface (plus van Hove singu- 
larities nearby) and does not include any information on 
the orbital content of the band at the Fermi level. Also 
numerical studies on the one-band Hubbard model us- 
ing the dynamical cluster approach^ and density matrix 
renormalization group calculations on t-J ladders^ have 
shown the trend that a larger t' leads to smaller Tc, in 
contradiction with the above mentioned findings. Only 
in t-t'-J, comprehensive Quantum Monte Carlo (QMC) 
studies by Spanu et ah-- showed a slight enhancement of 
pairing at optimal doping with nonzero t', however, to 
a smaller extent than in earlier variational QMC studies 
by Shih et alM Hence in the one-band model, both at 
weak and strong coupling, the theoretical expectations 
are inconclusive and certainly not fully consistent with 
the empirical trend. The simple question now is whether 
the theoretical picture is altered when orbital information 
of the multiband case is included. Therefore we revisit 
this problem in simple multiband models for the band 
structure of the cuprates. 

The change of the next-to-nearest-neighbor hopping 
parameter t' is caused by changes in the multiorbital 
electronic structure of the cuprates. In downfolded four- 
orbital models for the cupratesj^ii^ its increase is related 
to a lowering in energy of the so-called axial orbital to- 
ward the Cu 3c?a.2_y2-level. The axial orbital is basically 
a linear antibonding combination of the local Cu 4s and 
the surrounding oxygen 2pz states. It can be decreased 
by reducing the overlap between these two orbitals, which 
happens if the oxygens move further out of plane. In this 
way, crystal and electronic structures are correlated, and 
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the hope is to relate the structural differences between 
different cuprates to the differences in their supercon- 
ducting properties. 

The material trend pointed out by Pavarini et al^ 
also spurred exciting suggestions to produce band struc- 
tures and Fermi surfaces with even higher r parameters 
in LaNiOa/LaAlOa-heterostructures,^^ For such systems, 
taking over the trend from the cuprates would result in 
Tc's above lOOi^T. Therefore, a more detailed understand- 
ing of the relation between the low-lying electronic struc- 
ture and the superconducting transition temperature be- 
comes an important question in the field of tailored tran- 
sition metal oxide systems. 

The relation of finer differences in the multiband elec- 
tronic structure to the pairing strength for different 
cuprates has been addressed theoretically by at least 
two works. Kent et al~ studied downfolded dpp three- 
band Hubbard models using the dynamical cluster ap- 
proach with QMC impurity solver. The authors found 
a very strong sensitivity of the resulting Tc on choice 
of the downfolding technique and the localization of the 
Wannier functions. Small longer-ranged hopping had a 
marked influence on the results. We take this as an indi- 
cation that a direct parameter-free theoretical approach 
with a nonperturbative (cluster) many-body technique is 
still too challenging, and that a qualitative understanding 
of how the different building blocks of the model affect 
the resulting T^s would be very useful. 

A theoretical study that was successful in obtaining 
a significant difference in pairing strength between La- 
based and Hg-based cuprates in the right direction came 
from the group around Kuroki, Arita, and Aokij^ These 
authors considered two- and three-orbital models ob- 
tained using maximally localized Wannier orbitals to 
represent the low-lying density functional theory (DFT) 
band structure. The models were then treated by 
the fluctuation-exchange approximation (FLEX), which 
yields effective coupling strengths for d-wave pairing. 
Upon changing the oxygen height ho with respect to the 
Cu plane, the band structure was altered continuously 
from a situation corresponding to La2Cu04 to roughly 
that of HgBa2Cu04. The parameter in the two- or 
three-orbital models, that responded most to this struc- 
tural change, was the difference AE of the onsite kinetic 
energy between the two Cu ?>d^2_y2 {d^2_y2)- and Cu 
3^3x2 _^2(d22)-dominated Wannier states. If a third Cu 
4s(s) orbital was included into the model, its energy dif- 
ference to the dz2 was found to be roughly constant along 
this patch, while in the two-orbital model, this change 
was effectively absorbed into the model parameters. In 
both models, the coupling strength for c?-wave pairing 
was found to increase when ho or AE were increased, i.e., 
when the axial s level moved closer to the (i2.2_j,2-level 
from above, and the moved further down below the 
d2;2_j,2-level. Along this patch the Fermi surface became 
more rounded. This study shows, that differences in pair- 
ing strength can indeed be related to changes in two- or 
three-band model parameters. The present study was 



motivated by this work. Our goal was to acquire a clearer 
understanding why these model parameters changes ac- 
tually the pairing strength in a d-wave pairing situation, 
and if one can identify a simple mechanism behind the ob- 
served trends. This might be useful for prescribing other 
band structures that should have high transition tem- 
peratures. Below we will show that in the weak-coupling 
picture, there is actually a parameter window, where the 
detrimental effect of orbital admixture to the da;2_j,2-like 
conduction band is reduced more strongly than the other 
negative factor of Fermi surface rounding increases. This 
can lead to an increase in the pairing scale as function of 
the relevant parameters. 

This paper is organized as follows. First, we describe 
the underlying models and the functional renormalization 
group (fRG) method used for the study of these models. 
Then, we apply the fRG to the two-band model and find 
a first trend. This trend can then be understood more 
deeply by using the simplified two-patch model. Finally, 
we then use the three-band description that allows us 
to qualitatively reproduce the trends seen in the FLEX 
calculation^^ and to establish a simple picture. 



II. DESCRIPTION OF THE MODELS AND THE 
SCHEME 



In this paper, we first study a two-orbital model that 
can be thought to arise from one s-likc and one pla- 
nar c?^2_y2-likc orbital on the two-dimensional square 
latticed All other relevant orbitals, e.g., oxygen p states, 
should be considered as included in these effective or- 
bitals. The free Hamiltonian reads 



H 




es{k) + AEsd- IJ- v{k) 

v{k) ed{k) - fi 




with the nearest-neighbor hopping dispersions eg/dik) = 
—2ts/d{coskx + cos ky), the hybridization term v{k) = 
kx cos /cy) and the chemical potential /i. Here, 
is/d denotes the hoppings within the s or dx2-y2 orbital 
respectively, tsd is the hopping between the two orbitals, 
k is the momentum vector in the Brillouin zone and cr 
denotes the spin projection, c^^'''"^'' is the annihilator 

(creator) of a single-particle excitations in the s or dx2-y2 
orbital respectively. The symmetry of v{k) is due to the 
symmetry of the orbitals that hybridize. It is zero in the 
Brillouin zone diagonal and strongest at (tt, 0), (0,7r). 
The lattice constant has been set to unity. 

The matrix in orbital space can be diagonalized by 
a unitary transformation, leading to two bands. The 

new operators a-^^'' that annihilate and create the single- 
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particle excitations in band b are given by 



o 



where the sum in the first Hne runs over the (1^2 _y2 a 
the s orbital and the sum in the second line runs O" 
the two corresponding bands. The matrix element of 1 
transformation at a given momentum vector k is denol 
by Ubo{k). In the case of the three-band-model gi\ 
in Eq. ([7]), the transformation additionally includes 1 
contribution from the orbital. 

For the creation operators, we have to use the adjo 
equations. If the hybridization is not too strong, we c 
talk about a s-dominated and a dj.2_j,2 -dominated bai 
We now consider the situation where AEgd is rather la: 
and positive, and only the (i2,2_y2-dominated band ha 
Fermi surface. For chemical potential ^ = and la: 
AEsd, the Fermi surface is a perfectly nested square, 
while if we decrease AEsd, the Fermi surface gets rounded 
by the hybridization term. For the lower band around the 
Fermi level, we have the dispersion 

Ec(k) = 



\ 



es{k) + AEsd ~ ed{k) 



+ v^k). (3) 



The Fermi surface opens up at the van Hove points and 
becomes similar on shape as the Fermi surfaces observed 
in high- Tc cuprates. In Fig. [T]we show two examples. In 
these plots, the maximal admixture of the s orbital to the 
band with Fermi surface is only about 10%, but below 
we will see that already this has a measurable effect on 
the critical scale. 

Now, we include interaction in the form of intraorbital 
and interorbital repulsions, Ug/d = U and U' as well as a 
Hund's rule coupling Jh and pair hopping term Jp. We 
write the interaction Hamiltonian as 

U' 



(4) 



The indices i and o, o' = {s, d} denote the lattice sites 
and orbitals respectively and n° ^ — c°^„c°^. More gen- 
erally we can write this interaction in wave vector space 
as 
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^2,^3) 



xc°^^ c°^^ cZ^ c°-} . 

k^,(T k4,(T' k2,(7 ki,(T 



(5) 
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FIG. 1: Left: Fermi surface of the two-band model for large 
AEsd = 12eV; td = 0.45ey, ts = 0.5eV, tsd = 0.6eV (upper 
plot), and the orbital weights of the lower d^2 _y2-\\kB band 
through the Brillouin zone (lower plot). The vertical compo- 
nent of the vector arrows represents the d^2_y2 orbital ad- 
mixture, and the horizontal component denotes the s orbital 
admixture. The hybridization is strongest near the (vr, 0)- and 
(0, 7r)-points. The filling is fixed to (n) — 0.84 per site. Right: 
The same for smaller AEsd ~ deV . 



where k^ is fixed by momentum conservation on 
the lattice and JV is the number of unit cells. 



01,02,03,04 



(fci, ^2, fca) has to be chosen appropriately, but 
is wave vector independent for local interactions (and has 
otherwise some simple k variation). We can now write 
this interaction using the band operators a'S^"^ . This gives 



k ]^ , fc2 7 ^3 1 1 
01,02,0^,0^ 



bi,b2,b3,bi 



01 ,02,03,04 (fcl 1 ^2 I ^3) 

Moibi iki)uo.,h2 (^2)^0363 (fc3)w*,b^ (ki) 



^, bsi 64 1 fco bi 
fc3,cr A;4,ct' fc2,cr' fci,cr 



1 

2Af 



E 

k ]^ , fc2 , fc3 , cr , (T 



feat fe4t fe2 fei 
X a-r ' a-, ' a- a- 



(6) 



The factor Uoibi (A:i)uo2fc2 (^2)^*363 (^3X^64 (^4) contain- 
ing the orbital content is also referred to as "orbital 
makeup." 

The orbital content is the new aspect in comparison 
with one-band models. In the multiband case, already 
the bare interactions in band language are wave vec- 
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tor dependent. In context with unconventional super- 
conductivity, this raises interesting questions, in particu- 
lar, whether the bare interaction already contains an at- 
tractive component in some symmetry channel, or if the 
generation of such an attraction (that in one-band mod- 
els is usually accomplished by particle-hole-fluctuations) 
is somehow influenced by this extra structure. In the 
Sec. |V] on the two-patch model, we will give more de- 
tails which interaction processes are increased by that 
and which are reduced. We will see that under quite 
general conditions, the orbital content does not help d- 
wave superconductivity. 

Below we will also study a three-band model, which 
now also contains an orbital below the Fermi level, with 
the symmetry of a d^a-like orbital. For this model, we use 
parameters obtained with the Wien2Wannier schemeji^ 
We drop all hoppings in the third direction, in order to 
keep the model two dimensional. The free part of the 
Hamiltonian is then given by 



/4, 



k.(7 



with the diagonal terms 



es(fc) hsd{k) hsz{k) 
hsd{k) ed{k) hdzik) 
hszik) hdz{k) ez{k) 





(0,0) 



FIG. 2: Left: band structure of the three-band model with 
S-, d^2_y2- and d^2-\ike orbitals, with parameters correspond- 
ing to La2Cu04. Right: same band structure for larger 
AE = — Q but same Iq q — tg ^, roughly correspond- 
ing to HgBa2Cu04 according to Wien2Wannier.— Filling: 
(n) = 2.84. 



Again we use interactions of intraorbital, interorbital, 
Hund's rule and pair hopping type according to Eq. 
where now the sums o and o' go over all three orbitals o = 
{s, d, z}. In order to keep the analysis simple, we refrain 
from using orbital dependent interaction parameters. 



III. FUNCTIONAL RG TREATMENT FOR THE 
CONDUCTION BAND AND ADDITIONAL 
SECOND-ORDER CORRECTION 



and the hybridization terms 

hou02ik)= Yl C:Xe*^"=='=-+'"«'='') , (9) 

where the coefficients „^ and f^f^y > with o, oi , 02 
denoting the respective orbitals, are the parameters de- 
rived from the ab-initio calculation. The indices s, d, 
and z are used to label the s , d^2_y2 , and dz^ orbitals re- 
spectively, to shorten the notation. We include hoppings 
between sites on the two-dimensional square lattice with 
a maximum distance of five lattice constants in each lat- 
tice direction, that is the sum in Eqs. ^ and ^ runs 
over all integer values rnxjiriy from —5 to -1-5. 

The resulting band structures are plotted in Fig. [21 
Now we have one s-dominated band above the Fermi level 
and one -dominated band below the Fermi level. The 
band at the Fermi level is primarily of d^2_j,2-character. 

When applying the fRG technique described below, we 
only integrate out the dr^2_y2-\\ke conduction band with 
the Fermi surface, as described in Sec. IIIII We use a 
basis in which the quadratic part of the Hamiltonian is 
diagonal and thus have to transform the interacting part 
according to Eq. ([5]). Thus the influence from the ^^2 
and s orbital is then included in the perturbation of the 
dispersion relation and the transformation of the inter- 
action, i.e. the orbital mixing of Eq. ([B]). 



In this paper, the relevant situation is when the dx2-y2- 
like band has a Fermi surface and the other s-like and dz^- 
like bands are energetically separated from the Fermi sur- 
face. Then, if we are interested in the interaction effects 
at low temperatures and set up a perturbation expansion 
in the bare interactions, we can expect that virtual ex- 
citations into the bands away from the Fermi level only 
play a smaller role, as they generally lead to larger energy 
denominators in the corresponding diagrams. The more 
important point where the multiband character comes in 
may be the orbital content in Eq. ^ that makes the bare 
interactions projected onto the lower band wave vector 
dependent. In the following we will be concerned with 
the question how these extra factors in front of the in- 
teraction constants affect the critical scales for d-wave 
pairing. 

The next step is to sum the perturbation series in the 
interaction, i.e., to integrate out the band that forms the 
Fermi surface. This is done using the fRG iV-patch ap- 
proach, detailed, e.g., in Ref. [ij and ITsl This scheme 
takes into account the fermionic single-particle excita- 
tions step by step, with a decreasing infrared cutoff A. 

This leads to flow of the effective interactions for quasi- 
particles below this scale that contain the multiple scat- 
tering processes with higher-energy intermediate parti- 
cles. The present approach makes three simplifying ap- 
proximations discussed in previous works j^^^ First self- 
energy is neglected in the flow. We expect the self-energy 
to become important only when the effective interac- 
tions become large, and we stop the flow when this oc- 
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curs. Then, second, the one-particle-irreducible six-point 
vertex is ignored. Such a term is zero in the bare in- 
teraction, but will get generated. Its feedback on the 
normal two-particle or four-point vertex contains either 
two-loop diagrams (for which there is small phase space 
near the Fermi level) or corresponds to self-energy cor- 
rections on internal lines, which are ignored. Third, the 
frequency dependence of the interaction is ignored and 
all vertices are studied at zero external frequencies. The 
wavevector-dependence is treated in the iV-patch approx- 
imation which resolves the angular dependence around 
the Fermi surface by introducing N sectors with Fermi 
wave vectors ki labeled by indices ki with i = 1 . . . N, 
but drops the radial dependence. All these simplifica- 
tions are basically dictated by feasibility, and there is no 
rigorous statement about their validity. However, there 
is also no clearly documented failure of these approxima- 
tions in reasonable context. Therefore the critical scales 
obtained this way should serve the more qualitative pur- 
poses of this paper. 

In this approximation, the flow equations for the cou- 
pling function V^^(fci, A:2, ^3) depending two incoming 
wave vectors ki/2 and one outgoing k^ read 



V^{ki,k2,k3) = T^p + T^H.d + rpH.cr , (W) 



dX 



with the particle-particle channel 

T^pp(fcl,fc2,/C3) 

^-TY,y\ki,k2,k)-L\k,qpp)V\k,qpp,k3), 

k,ikQ 

(11) 

the direct particle-hole channel 

= -TJ2 [-^y\ki,k, k3)L\k, qPH,d)V\qpH,d, k) 

k,ikQ 

+ V\k, kuk)L^{k, qPH.d)V^{qPH,d, k2,k) 

+ V\kuk,k3)L\k,qpH,d)V\k2,qPH,d,k)] , (12) 

and the crossed particle-hole channel 

PH,criki,k2, fcs) 
= _ r ^ V\k, k2, h)L\k,qpH.,,r)V\h,qPH,cr, k) , 

(13) 

where k — {ikQ,k) are wave vector and Matsubara fre- 
quency of the first line and qpp = ~k + ki + k2, qpH,d = 
k + ki — ks, qpH,cr = k + k2 — k3 are the wavevectors of the 
second loop line. T is the temperature. All incoming and 
outgoing frequencies are set to zero and the frequency of 
the second line is fixed by frequency conservation to be 



—iko in the PP diagram and iko in the PH diagrams. 
Also the fourth momentum in the interaction vertex is 
fixed by conservation, and the spin convention for the 
coupling function A;2, fcs) is that the first incom- 

ing line fci and the first outgoing line ^3 have the same 
spin projection a, while the second incoming k2 and the 
second outgoing have spin a' . The internal loop is given 

by 



L\k,k') = -[G\k)G\k')] 



(14) 



where in our approximation self-energy corrections are 
neglected, i.e., the full propagator is identical to the free 
propagator. The fRG cutoff is introduced by multiplying 
the Green's function by a function that cuts out the low- 
energy modes. 



G\k) = 



C^[E{k)] 
iko - E{k) 



(15) 



In the numerical implementation of the flow scheme, we 
use a slightly smoothened step function, C''^[E{k)] « 
O ^|i?(fc)| — . The initial condition for the flow at the 

scale Ao of the band width of the conduction band is 
given by the orbital-content dressed coupling function in 
Eq. (0), i.e., by 



V^"{ki,k2,k3)^Vccccikl,k2,k3). 



(16) 



with band index c for the c?2.2_j^2-dominated conduction 
band. 

The integration of the above fRG equation takes into 
account all one-loop corrections with both internal lines 
in the low-energy window and thereby reconstructs the 
most important parts from the perturbation expansion 
of the band near the Fermi level. Usually, the fRG flow 
leads to strong coupling, i.e., at some low critical energy 
scale Ac, at least one family of couplings V^{ki,k2,k3) 
grows very large, and the perturbative flow without self- 
energy feedback has to be stopped. The different types of 
divergences for the single-band Hubbard model as func- 
tion of the system parameters have been analyzed and 
classified in a number of previous publications (see, e.g., 
Refs. 17 and 18). In the present case the fRG flow al- 
most exclusively favors a d-wave pairing instability that 
is understood to be driven by antiferromagnetic spin fluc- 
tuations. In this paper we do not discuss any other fea- 
tures of this c?-wave instability except its critical scale 
Ac. This scale can be understood as a measure or up- 
per estimate for the transition temperature Tc into the 
d-wave superconducting state. It takes an analogous role 
as the quantity fiijjije^^''°^ in a simple BCS problem with 
Debye frequency wd, density of states po and attractive 
interaction g. 

So far, the bands away from the Fermi surface have 
only entered through the orbital content. Let us next 
discuss the first step to include virtual processes into 
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these bands, i.e. inter-band transition. If we write a per- 
turbation series for the full two- or three-orbital model, 
and compare the various one-loop terms entering there to 
what we get in the fRG for the conduction band, we find 
that the next important terms that are not included in 
the fRG treatment are diagrams with at least one inter- 
nal line in one of the bands away from the Fermi surface. 
While these diagrams will not become singular as their 
energy denominator will never go to zero, they might 
induce additional wave vector dependencies to the in- 
teractions, that might be measurable for the unconven- 
tional pairing studied here. Actually, the constrained- 
RPA formalism sums up a part of these contributions 
to infinite order, which results in an additional screen- 
ing of the Coulomb interaction.—-^ Including these ef- 
fects more systematically into the fRG for a single band 
is a separate issue that we do not touch here. Instead, 
we will study the impact of an additional second-order 
correction due to these virtual processes with one and 
also two internal lines away from the conduction band, 
whereas the contribution of the latter is confirmed to be 
comparably small. This correction, which we will refer 
to as high- energy second- order correction, is added to the 
initial condition of the fiow in the conduction band. It 
is given by all one-loop diagrams (particle-particle and 
particle-hole) with bare interactions at the vertices and 
at least one internal line in a band above or below the 
conduction band. We will see that it has a definite impact 
for smaller energy gaps between the bands. Due to the 
energetic separation of the bands away from the Fermi 
surface, we expect that this second-order treatment is al- 
ready quite good, and that the higher-order corrections 
neglected here only have a small effect. 



IV. RESULTS FOR THE TWO-BAND MODEL 

Let us now apply the A^-patch fRG scheme to the lower 
band of the two-band model described in Sec. II. We 
therefore fix the hopping parameters of the bands to be 
td — 0A5eV, ts = 0.5eV, and the hybridization to tgd — 
O.GeV. The tuning parameter is the energy separation 
between the band centers, AEsd- 

In the plots of Fig. [3] we additionally choose two 
sets of interaction parameters according to the sum rule 
U -U' = 2Jh = 2Jp^ The density is fixed to 16% hole 
doping of the lower band, i.e. (n) = 0.84 per site. In the 
left plot, the critical scale for the d-wave pairing insta- 
bility, Ac, in dependence of the energy difference AEsd 
of the two orbitals is shown. The direction of the hori- 
zontal axis is inverted so that the plots can be compared 
better with the ones for the three-band model in Sec. 
IVIl We can see that moving the two bands closer to each 
other, i.e., reducing AEgd, decreases the critical scale Ac. 
This behavior is robust with respect to the interaction 
strength. We also plot Ac vs AEsd when the orbital con- 
tent is ignored and the bare interaction is simply U for all 
wavevector combinations. In this case, Ac is higher. This 
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FIG. 3: Critical scale (left) and superconducting suscepti- 
bility at a fixed scale O.OQeV (right) vs AEsd, computed 
with 64 patches on the Fermi surface. Red triangles: with 
orbital mixing, green diamonds: orbital mixing cut out by 
hand, U = l.SeV, U' = l.UeV, Jh = Jp = 0.l8eV. Blue 
squares: with larger interactions U — 2.2eV, U' — 1.76eV, 
Jh = Jp = Q.22eV , with orbital mixing. Parameters: 
td = 0.45el/, ts = 0.5el/, tsd = 0.6el/, in) = 0.84. 



allows us to distinguish the effect of the orbital content 
from that of the Fermi surface shape, which also changes 
as a function of the energy separation AEgd- Here, we 
only calculated data for the parameter set with smaller 
interaction, as the one-loop approximation is not justi- 
fied, if the interaction is to large. With included orbital 
mixing however the average of the coupling function at 
the beginning of the flow is below half of the bandwidth. 
Summarizing our findings, we can state that the expected 
decrease of Ac with the Fermi surface becoming rounder 
when AEsd becomes smaller is not compensated by the 
orbital content or stronger hybridization in the d^^-y^- 
band. Rather, this latter effect leads to an additional 
reduction of the pairing scale. 

We also plot in Fig. [3] the effective strength of the 
d-wave superconducting channel at a small fixed A > 
Xc{AEsd), which is defined as 



Xsc 



±J2v\k,~k,p)f{k)f{p). 



(17) 



Here, the sums run over all A^ patches of our patching 
scheme and f{k) and f(p) are d-wave form factors and 
are given by f{k) — cos{kx) — cos{ky), and V\ (jt, —k,fj 

is the running interaction at the small fixed scale. This 
measure of the strength of the superconducting channel 
gives a qualitatively similar picture. Again, both factors, 
Fermi surface deformation and orbital makeup, reduce 
the tendency toward d-wave pairing. 

In this sense, a well separated (ia;2_j,2-like band appears 
to be the optimal situation for a high critical scale or 
high Tc- The two-band model with the s orbital above 
the d^2_y2 orbital can not give an explanation why the 
cuprates with rounder Fermi surface have higher Tc, as 
indicated by the DFT trends. Our data clearly show that 
the perturbation of the dispersion as well as the orbital 
content reduce the critical scale. 

In order to understand this finding more directly, and 
to see how general it is, we will now turn to a simplified 
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description in terms of the two-patch model. 



V. SIMPLIFIED PICTURE: TWO-PATCH 
MODEL 

Let us consider again the situation where only the 
lower d2.2_j,2-like band has a Fermi surface, and let this 
Fermi surface be near the van Hove points at A — (tt, 0) 
and B = (0,7r). Then a common approximation is the 
so-called two-patch modelj^^ where only the fermionic de- 
grees of freedom in small circles around these van Hove 
points are kept, and the interactions between these re- 
gions are approximated by four constants gi to 34. These 
are defined by 

gi=V{A,B,B,A)^ViB,A,A,B), (18) 

g2^V{A,B,A,B)^V{B,A,B,A), (19) 

g3^V{A,A,B,B) = ViB,B,A,A), (20) 

g4^ViA,A,A,A) = ViB,B,B,B), (21) 

where the notation, is again, that the first two entries 
ki and fc2 of V{ki, k2,ks, ^4) are the incoming wave vec- 
tors, and that the first incoming particle with ki and the 
first outgoing with fca have the same spin projection a, 
while the second incoming ^2 and the second outgoing 
^4 have a'. Such a modeling has been used to explore 
the basic phase diagram of the two-dimensional single- 
band t-t' Hubbard model near van Hove filling. There, 
the initial value for the gi from the bare interaction is 
just U. Recently, a very similar model has been used for 
the study of the iron arsenide superconductors.^^ 

The fRG treatment of the two-patch model for the 
Hubbard model, where only the singular one-loop contri- 
butions of the particle-particle loop with zero incoming 
momentum and the particle-hole loop with momentum 
transfer (tt, tt) are kept, leads to the fiow equations 



91 

92 
93 
94 



2digi(52-3i), (22) 

di{9i+9l), (23) 

-2^05334 + 2(iig3(252 - 9i) , (24) 

-M9i+9l), (25) 



with do,di > as scale derivatives of the particle-particle 
and particle-hole loop, respectively. Possible instabili- 
ties in the d-wave pairing, antiferromagnetic spin-density 
wave (SDW), and charge-density wave channel are then 
associated with the unbounded growth of the combina- 
tions 33-54, 52 + 53 and ~2gi +52-53, respectively. 
In the one-band Hubbard model for do = di , correspond- 
ing to perfect nesting, the coupling constants diverge as 
52,3 00, 54 —00, gi diverges more slowly. As long 
as 54 is positive, the two terms on the RHS of Eq. (p4| 
have opposite sign and give competing contributions, but 
eventually 54 becomes negative and both terms in Eq. 
(1^^ drive the flow to strong coupling. Then the d-wave 



pairing g^ — g^ diverges, but usually, for nonzero U, the 
SDW combination 52 +53 has grown much larger. Hence 
the instability for di = do is of multichannel type, with 
the SDW tendencies being strongest. For di < do, cor- 
responding to nonzero t', the SDW channel is less domi- 
nating and one gets a regime where the d-wave coupling 
53 — 54 diverges more strongly. This is what we call the 
d-wave pairing regime. A more detailed analysis is given 
in Ref. M- 

For the single-band model, the comparison with the 
fRG calculations taking into account the full wave vector 
dependence of the interactions around the Fermi surface 
showed that most trends can already be inferred from 
the two-patch model, so that the latter should serve as a 
good starting point to study the main effects. 

We will now use the two-patch model to assess the 
impact of the orbital content on the type and energy scale 
of the instability of the weakly-coupled state due to the 
interactions. At the van Hove points, the transformation 
that diagonalizes the hopping Hamiltonian of the two- 
band model can be expressed as 



UcdiA) = COS0, UcdiB) = cos(f>, 
UcsiA) = sin0, Ucs{B) = - sint^ 



(26) 
(27) 

where Ucs and Ud are the orbital weights of s and dx2_y2 
orbitals in the conduction band c, which is the one with 
the Fermi surface near the van Hove points. The angle 
(/) is a measure for the strength of the orbital mixing 
and can be understood as the angle between the vertical 
and the arrow at the van Hove points in Fig. [TJ In the 
case oi <j) = 0, the mixing is zero and we recover the 
pure one band model without mixing effects. The sign 
change of the s admixture in Ucs{B) compared to Ucs{A) 
comes from the different in-plane symmetry of the dx2-y2 
orbital. It is crucial for the following. Note that this 
sign change would also occur for any other orbital that 
is admixed to the dx2_y2 orbital that transforms trivially 
under 90 degree rotation in the plane, i.e. also for d^^ or 
Pz orbitals. Hence, qualitatively, the results hold more 
generally for a wider class of admixed orbitals. 

Next we consider here only local intra- and intcror- 
bital repulsion U and U' < U to keep it simple. Hund's 
couplings and pair hopping terms in a realistic range are 
smaller than these parameters, and should not lead to 
qualitative differences in the results. For the intraor- 
bital interaction, the orbital content gives a factor cos^ (j), 
while for the interorbital interaction U' the two incoming 
and the two outgoing particles are in different orbitals, 
contributing ±2 sin^ (j) cos^ (f>. The sign now depends on 
whether the particle in the s orbital remains in the same 
patch, say A, or whether it gets scattered to the other 
patch, say from A to B. In total, we obtain for the initial 
values of the four coupling constants: 



51,3 
52,4 



U 
U 



-~2U'cos^(j)sm^(j), (28) 
+ 2U'cos^(j)sm^(j). (29) 



We see that the orbital content (/> 7^ in this d3,2_j,2-case 
quite generally suppresses the bare interaction strength. 
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FIG. 4: Critical scale of the two-patch model as a function 
of the mixing parameter (j) for [7 = 1 and different U' . For 
most parameters, the instability is of antiferromagnetic spin- 
density wave or d-wave pairing type, and the critical scale 
decreases with larger admixture <j). Only for U' > U and 
(/!> large enough, the instability turns into a charge-density 
wave instability, and the critical scale is increased by further 
increasing (p (black dashed dotted curve). 

For a nonzero value of U' > 0, this suppression is stronger 
for (7i_3 and weaker for 172, 4- Generally, a larger 172 and 33 
would lead to a higher critical scale for d-wave pairing, as 
these couplings drive the antiferromagnetic fluctuations 
that form the pairing glue. More repulsive gi and (74 leads 
to a smaller critical scale. For the desired d-wave pairing 
instability, the relative enhancement of 54 with respect 
to 53 for ^ is exactly the wrong way. It corresponds 
to an additional repulsion in the d-wave pairing channel 
and reduces the critical scale. 

We have performed a numerical analysis of the flow 
equations ([22)) -p5 |) . in order to investigate the influence 
of the mixing on the critical scale. Fig. U shows the 
critical scale obtained for the two-patch model for differ- 
ent [/'. It turns out that the critical scale is suppressed 
by mixing as long as U' does not get too large. Only for 
U' > U, we encounter the charge-density wave instability 
and here the s admixture — ^ 7r/2 can actually increase 
the critical scale in a certain range. Here, the initial value 
of gs becomes negative and thus 53 flows to —00, as the 
second term of the right hand side of Eq. (|24p is nega- 
tive and the first term eventually becomes negative when 
54 < 0. In this numerical treatment, we assumed that 
the relevant loops do and di have the same value during 
the flow, which should be an acceptable assumption for 
these qualitative conclusions. However we checked that 
the qualitative results do not depend on the individual 
behavior of the loops during the flow. 

Thus we conclude that already the simple two-patch 
model can give us a simple explanation how orbital mix- 
ing effects lower the critical temperature of a d-wave su- 
perconducting instability. 

VI. RESULTS FOR THE THREE-BAND MODEL 

So far, we have seen that admixing a single orbital 
to the da.2_„2 -derived band leads to a reduction of the 



critical scale for d-wave pairing. This causes problems 
to understand the mentioned apparent material trend in 
the high- Tc cuprates that lowering an axial orbital in en- 
ergy from above results in higher critical temperatures. 
Here, we will show that this trend can be recovered if 
additional orbitals are considered. To this end, we now 
consider a two-dimensional three-band model, which in- 
cludes the dj,2_j^2, dz2, and s orbitals. The Hamiltonian 
of this model is given in Eq. ([7]). 

To find out whether the mentioned material trend can 
be explained by mixing effects, we study the critical scale 
for d-wave pairing with different values for the onsite ki- 
netic energy of the (1^2 and s orbitals. We follow the 
path through the parameter space proposed by Sakak- 
ibara et al^ More precisely, we increase AE = ~ ^0 
by hand but leaving the difference between d^2 and s 
orbital constant. This is motivated by the observation 
that in the cuprate HgBa2Cu04 with Tc ^ 90K both the 
and s orbital are lowered with respect to the d^^-y-i 
orbital compared to La2Cu04 with lower Tc ^ 30K. 
HgBa2Cu04 has a higher Tc and a more rounded Fermi 
surface than La2Cu04. As the change in AE appears 
to be the most striking effect, we do not consider a vari- 
ation of other hopping parameters to keep our analysis 
simple. Note that a larger AE = ~ ^0 ^^^^ corre- 
sponds to a smaller AEgd of Sec. llVi which is written as 
AEsd ~ ^0 ^0 

g in terms of the three-band model in 
Eqs. ([T])-©. All calculations presented in this section 
are done at 16% hole doping, that is, (n) — 2.84 per site. 

The left plot of Fig. [S] shows the Fermi surface of our 
model for the two situations. While the s orbital lowers 
the energy of the van Hove point, the dj;2 orbital has the 
opposite effect because it lies below the dx2_y2 orbital. 
Thus the combined effect of the s and dz2 orbitals can 
lead to a roughly squarelike Fermi surface, as it would 
be without mixing at all. The plot on the right shows 
the relative orbital admixture at the van Hove point for 
different AE. \uco{tt, 0)p is the squared absolute value of 
the coefficients in the orbital-band transformation given 
in Eq. where the index c as above denotes the da,2_y2- 
dominated conduction band. This quantity serves as a 
measure for the orbital admixture. The contributions 
from the three orbitals add up to one due to the unitarity 
of the transformation. With larger AE the Fermi surface 
becomes more rounded, due to a larger perturbation from 
the s orbital and a smaller influence from the dz2 orbital. 
The total mixing, i.e., the non-da;2_j,2-content, is weak- 
ened with larger AE. This mainly comes from the de- 
creasing influence of the dz2 orbital when it moves down 
in energy. On the other hand, the s orbital gets closer to 
the dx2-y2 orbital with larger AE and thus increases the 
mixing. This effect, however, is much smaller due to the 
larger gap between dx2_y2- and s orbitals, and because 
the s band is rather wide. We conclude that the lower- 
ing of AE leads to two competing effects on the critical 
scale. First, the Fermi surface gets more rounded, which 
according to the results in one-band models decreases the 
critical scale. Second, and this is now different from the 



FIG. 5: Left: shape of the Fermi surface for parameters 
corresponding to La2Cu04 (blue sohd hne) and for AE ~ 
2eV (red dashed hne), corresponding to HgBa2Cu04, (n) — 
2.84. Right: admixture of the three orbitals to the d^2_y2- 
dominated band with the Fermi surface, the orbitals d^-2_y2, 
d^2 and s are shown as a blue solid, green dashed, and red 
dotted lines, respectively. 



0.03r 6 




A E/eV A E/eV 



FIG. 6; fRG data for the three-band model with interaction 
parameters U = 2.2eV , U' = 1.76eV, Jh = Jp = 0.22eV 
computed with 64 patches, (n) = 2.84. Left plot: criti- 
cal scale Ac vs AE = fo,o — to,o- Right plot: average d- 
wave pairing at a fixed scale A = 0.027eV^ vs AE; without 
high-energy second-order correction (blue squares), with high- 
energy second-order correction (red triangles), explained in 
Sec. III. 



two-orbital case, the orbital mixing is reduced, and the 
dj:2_y2 character increases. As we have seen in the pre- 
vious section, this favors an enhancement of the critical 
scale. 

In Fig. [HI the critical scale, obtained from the fRG 
calculation, is shown. Here, the interaction parameters 
are chosen to he U — 2.2eV, Jh ^ Up ~ O.IU, and 
U' — 0.8U. We use 64 patches to discretize the Brillouin 
zone according to the patching scheme described above. 
With the interaction parameters chosen here the average 
coupling constant at the beginning of the flow is below 
half of the bandwidth. Both approximations, with and 
without high-energy second order correction, show the 
same qualitative trend. The two competing effects now 
induce a maximum in the critical scale and xsc a-t some 
AE*. For larger AE > AE* the reduced orbital mixing 
of the vertices cannot compensate the change of the Fermi 
surface. To separate the effect of the orbital mixing on 
the interaction and the dispersion, we run the same fRG 
flow, but now starting with a simple onsite repulsion re- 
gardless of the orbital mixing. In other words, we cut out 
the mixing effect on the vertices by hand. As expected 



FIG. 7: fRG data for the three-band model with interaction 
parameters U = 2.0eV, U' = 1.6eV, Jh ^ Jp = 0.2eV 
computed with 64 patches, (n) = 2.84. Left plot: criti- 
cal scale Ac vs AE = fg.o ~ ^o,o- Right plot: average d- 
wave pairing at a fixed scale A — 0.015eV^ vs AE; without 
high-energy second-order correction (blue squares), with high- 
energy second-order correction (red triangles), explained in 
Sec. III. The steps in the data are due to the discretization 
of the wave vector dependence of the coupling function. 



from the results of the two-patch model the critical scale 
without the orbital content is indeed drastically increased 
and we do not get an enhancement of the critical scale 
with increasing AE, because now only the Fermi surface 
effect of the shift comes into play. 

Additionally we run the fRG flow again but now shift 
only the s or the d^^ orbital respectively to see, if the 
interpretation suggested from the two-patch-model holds. 
As expected, if we only lower the s orbital, the critical 
scale decreases, since the mixing of the orbital is not 
reduced and the positive effect on the critical scale is 
not present. If we on the other hand lower only the d^^ 
orbital but keep the s orbital constant, the qualitative 
picture from Fig. [5] is recovered. 

The value AE* with highest critical scale depends on 
the strength of the interactions. For smaller bare values, 
the Fermi surface shape is more important, and AE* 
shifts to smaller values and even disappears. In Fig. [71 
we again show critical scale and d-wave pairing strength 
vs AE in Fig. [SI but now for reduced onsite interac- 
tion U = 2eV , while the other interaction parameters 
are reduced accordingly so that the ratio to the onsite 
interaction remains constant. Now without second-order 
correction of the additional orbitals the critical scale de- 
creases over the entire range. Our data quite generally 
suggests that at least without high-energy corrections the 
effect of an enhancement of the critical scale with larger 
AE is less pronounced, when the critical scale is lower, 
that is for smaller interaction strength or, e.g., larger 
J/f , which usually suppresses the flow to strong coupling. 
Turning this observation around, we can argue that the 
enhancement effect might even be larger if we used real- 
istic interaction values for the cuprates. We refrain from 
running the fRG flow with such large interaction param- 
eters, as the one-loop approximation of the fRG equation 
is then no longer justified. 

It can be seen that the inclusion of the high-energy 
vertices within second-order perturbation can have a 



10 



profound effect on the critical scale. Especially for 
smaller interaction strength, the instability is strongly 
suppressed at small AE, so that in contrast to the case 
without high-energy correction even at smaller interac- 
tion strength a maximum in the critical scale can be 
found. It seems worthwhile to study the renormalization 
of the vertices due to higher energy modes in more detail. 
This may be an additional contribution to the observed 
material trends, as the curves with included high-energy 
corrections exhibit a stronger increase in the critical scale 
than the ones without. 

We can compare our findings and our idea, why this 
Tc-trend occurs, with the FLEX two-orbital and three- 
orbital calculations by Sakakibara et al^ Although their 
Eliashberg eigenvalues cannot be compared directly to 
the fRG critical scales, it seems that we can qualitatively 
reproduce their results. Clearly, at least for their two- 
band calculation, the orbital admixture to the conduc- 
tion band is certainly lower for the higher Tc's with the 
rounder surfaces, supporting our explanation. Quantita- 
tively, our results are sensible to the interaction strength 
and a direct comparison would require data at larger on- 
site interaction of ieV as used for the FLEX calculation. 
As the band width of the conduction band is only ~ 4eF 
we refrain from running the fRG for such large interac- 
tions. Also we expect that self-energy effects, which were 
included by Sakakibara et al. but not in our work, will 
probably still cause a noticeable quantitative difference. 

Within this downfolded three-band electronic struc- 
ture, HgBa2Cu04 with Tc of 90K would roughly corre- 
spond to AE = 2eV, close to where the maximal criti- 
cal scale occurs in Fig. [B] While this is promising, we 
note that the band structure change from La2Cu04 to 
HgBa2Cu04 may be more complex, and taking into ac- 
count these additional changes may affect the observed 
enhancement as well. In any case, on a quantitative level 
and possibly different from Ref. our weak-coupling 
fRG studies cannot be expected to explain the experi- 
mental trend, because the real cuprates are more strongly 
interacting. We propose that strong-coupling methods 
should be used to find out whether the picture drawn 
here works the same (and quantitatively better) at strong 
coupling. Nevertheless, we have identified AE' as an im- 
portant tuning parameter. 

The model with the La2Cu04-like band structure pa- 
rameters has a nested Fermi surface. Consequently, aside 
from the superconducting instability, we observe a strong 
tendency toward long-range antiferromagnetic (AFM) or- 
der. With larger AE, the nesting of the Fermi surface 
becomes poorer, so that the AFM tendency is more and 
more suppressed. To have a closer look on the two com- 
peting instabilities, we compare the effective strengths 
of the respective channels. The d-wave superconducting 
channel can be measured by xsc given in Eq. ()17p . Sim- 
ilarly, the AFM channel is 

XAFM = ^J2^^ ^2 « i^, 0), fca « (0, ^)) , (30) 

k 




log(i) log(i) 

FIG. 8: Flow of the antiferromagnetic (green dashed line) 
and superconducting (blue solid line) interaction strengths 
for the La2Cu04-like band structure (left plot) and for the 
HgBa2Cu04-like band structure with AE ~ 2eV (right plot) 
for the same parameters as in Fig. [S] In the first case, both 
channels compete, while in the second case, the d-wave pairing 
channel clearly dominates at low scales. 

where k2 and ^3 are chosen so that they lie close to the 
two van Hove points and the corresponding momentum 
transfer is k2 — ~ (t, tt). The corresponding order- 
ing susceptibilities diverge if these averaged strengths of 
the channels diverge. Thus the defined quantities can be 
regarded as a measure for the strength of the respective 
instabilities. The left plots of Figs. [S] and [7] show that 
the superconducting channel has, as expected, a similar 
behavior as the critical scale as function of AE. However 
atU — 2eV a qualitative difference arises in the effective 
strength of the superconducting channel in comparison 
with the critical scale (see Fig. [7]) . The effective strength 
still exhibits a maximum, although it is not as distinct 
as for larger interactions. 

We also compare the relative strength of the channels 
for different AE in Fig. |51 For small AE in the left 
plot the AFM tendency is larger than the d-wave pairing 
strength, and both averages grow with similar exponent. 
As AE becomes larger, the AFM tendency is weakened 
due to poorer nesting of the Fermi surface. Then the 
superconducting channel is leading at low enough scales. 

Usually, the critical scale Ac is regarded as an upper 
estimate for actual transition temperatures into ordered 
states. Regarding the question what the true supercon- 
ducting transition temperatures might be in the two situ- 
ations compared here, it is now very likely that in the case 
of strong competition between two channels, the transi- 
tion temperature will actually be significantly reduced 
compared to Ac, while for the cleaner pairing instabil- 
ity in the non- nested situation, Tc might be closer to Ac. 
We expect that an inclusion of the self-energy, which is 
not considered in this work, would capture this effect. 
Thus the true superconducting transition temperatures 
at small AE would be even smaller, and consequently 
increasing AE should increase the superconducting chan- 
nel even more than in our calculation. 

Summarizing these observations we state that the mul- 
tiorbital model with orbital energy variations in accor- 
dance with actual material differences is able to repro- 
duce to some degree the phenomenological tendency that 
critical temperatures can grow although the Fermi sur- 
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face gets more rounded. This trend cannot be understood 
within the framework of one-band models. At least on 
a qualitative level, we have shown that the orbital mix- 
ing provides a mechanism for a enhancement at weak 
to moderate interactions. If this mechanism still works 
at the realistic interaction strength, it could, at least in 
part, be responsible for observed material trend for the 
Tc's of the high Tc cuprates. 

VII. CONCLUSIONS 

We have studied two-orbital and three-orbital models 
on the two-dimensional square lattice that mimic the low- 
lying electronic structure of the high- Tc cuprates. Using 
the fRG, we have computed the critical scale for d-wave 
pairing, which can be used as an estimate for the critical 
temperature for d-wave superconductivity. In one-band 
models with simple onsite repulsion, this energy scale is 
mainly dominated by the Fermi surface shape and de- 
creases when the Fermi surface gets more rounded. The 
material trend for the real high-Tc cuprates, or more pre- 
cisely, the combination of experimental Tc and electronic 
structure calculations for a series of cuprates, seems to 
contradict this trend, as materials with more rounded 
Fermi surfaces have higher experimental Tc's. Our goal 
was to see if orbital admixture to the conduction band 
reverts the shape-related Tc variation and allows one to 
understand this material trend. 

From our studies, we can draw two conclusions. First 
of all, for spin-fluctuation-induced d-wave pairing on the 
square lattice at a higher critical scale, it is best to have 
an energetically well-separated dx-^-y^-Yike band. We 
have shown that the admixture of bands with symmetry 
different from d^^-y^ from above and below in energy 
typically reduces the critical scale for pairing compared 
to the situation without admixture. Reducing the con- 
sideration to the two-patch model allowed us to relate 
the orbital admixture to an additional repulsion in the 
d-wave channel, which explains the reduction of the d- 
wave pairing tendencies. In this setup it can also clearly 
seen that admixing orbitals of s, , or, possibly. Pz char- 
acter have a similar negative effect, as the sign structure 
of the admixture is the same. 

With the knowledge that a single nearby band disturbs 
the pairing, we then analyzed situations with more than 
two bands. Here, a change of the orbital energies with 
one empty band moving closer toward the Fermi level can 
still result in a relative enhancement of the superconduct- 
ing instability scale, consistent with the actual material 
trend. We have shown that in three-band models for the 
cuprates, the approach of the wider s band can have a 
smaller negative effect on Tc than the positive effect of the 
d^-^ band moving further down to lower energies simulta- 
neously. Hence we have identified a possible path how the 
(i-wave Tc, generated by the spin-fluctuation mechanism, 



can be increased as a function of the energetic separation 
of orbitals near the Fermi level, although the Fermi sur- 
face of the (ij.2_,y2 -dominated band gets more rounded. 
The main reason for this effect is that the Tc-raising re- 
duction in the orbital admixture to the conduction band 
overcompensates the decrease in Tc due to the rounder 
Fermi surface. 

We have seen that choosing DFT-derived band- 
structure parameters we can qualitatively reproduce ex- 
perimental trends in Tc differences between La2Cu04 and 
HgBa2Cu04. However our calculations are only valid in 
the weak to moderate coupling regime, and if the inter- 
action parameters are chosen to small the expected be- 
havior is not rigorously reproduced. On the other hand, 
we have shown that the Tc enhancement works better 
for larger interactions. The real cuprates are of course 
quite strongly correlated, beyond the interaction range 
where we can apply our method. There, the Tc enhance- 
ment might be even more drastic, but based on our weak- 
coupling study, this is far from proven. Furthermore, 
the Tc increase may be even stronger than indicated by 
our numbers already for weaker interactions, as the less 
rounded Fermi surfaces have a stronger channel compe- 
tition, which should result in a further lowering of the 
critical temperature for d-wave pairing. Both these ob- 
servations suggest that the mechanism considered in this 
work might be more effective in the actual materials. It 
should be interesting to look for the same tendency, e.g., 
with cluster-DMFT techniques. It should also be inter- 
esting to quantify the degree of admixture to the dx2_y2 
band through the series of high-Tc cuprates that obeys 
the proclaimed relation between Tc and the electronic 
structure, and to see whether the Tc increase is corre- 
lated with this. 

In this work, we have focused on one single aspect that 
may distinguish different cuprate materials, expressed by 
the orbital energies in three-orbital effective Hamiltoni- 
ans. This allowed to obtain some understanding of how 
the pairing tendencies change. Very likely, there are other 
model parameters that exert additional influence on the 
energy scale of superconducting pairing. For example, 
depending on the hybridization with the surrounding or- 
bitals and the spread of the respective Wannier functions, 
the interaction parameters might vary. It appears to be 
an interesting topic to study these additional effects again 
in isolated form, and to understand their importance. 
Then it might be possible to compose a combined pic- 
ture that might ultimately be used to guide the search 
for higher transition temperatures. 
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